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Abstract 

The continuous evolution of a wide variety of systems, including 
continous-time Markov chains and linear hybrid automata, can be 
described in terms of linear differential equations. In this paper we 
study the decision problem of whether the solution x{t) of a sys¬ 
tem of linear differential equations dx/dt = Ax reaches a target 
halfspace infinitely often. This recurrent reachability problem can 
equivalently be formulated as the following Infinite Zeros Problem: 
does a real-valued function / : R>o —>■ R satisfying a given linear 
differential equation have infinitely many zeros? Our main decid¬ 
ability result is that if the differential equation has order at most 
7, then the Infinite Zeros Problem is decidable. On the other hand, 
we show that a decision procedure for the Infinite Zeros Problem 
at order 9 (and above) would entail a major breakthrough in Dio- 
phantine Approximation, specifically an algorithm for computing 
the Lagrange constants of arbitrary real algebraic numbers to arbi¬ 
trary precision. 

Categories and Subject Descriptors F.2.m [Analysis of Algo¬ 
rithms and Problem Complexity }: Miscellaneous 

Keywords linear dynamical systems, reachability, differential 
equations, Diophantine Approximation, Skolem problem 

1. Introduction 

A simple type of continuous-time system is one that satisfies a 
linear differential equation ^ = Ax, where A is an n x n matrix 
of rational numbers and x{t) G M" gives the system state at time t. 
In particular, such differential equations describe the state evolution 
of finite-state continuous-time Markov chains (via the so-called rate 
equation) and the continuous evolution of linear hybrid automata. 

A fundamental reachability question in this context is whether 
x{f) infinitely often reaches a target hyperplane v"‘"x = 0, where 
V G R" is the normal vector. Such a hyperplane could represent a 
transition guard in a hybrid automaton or a linear constraint on the 
state probability distribution of a continuous-time Markov chain 
(e.g., that the probability to be in a given state is exactly one half). 

A solution x(t) to the differential equation ^ = Ax ad¬ 
mits a matrix-exponential representation x{f) = e^*a:(0). The 
problem of whether x{f) reaches the hyperplane = 0 in- 
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finitely often then reduces to whether v^x{t) = 0 for infinitely 
many values of f > 0. Now the function / : R>o —>■ R defined 
by /(^) = v^x{t) can be written as an exponential-polynomial 
f{t) — , where the aj and the coefficients of 

each polynomial Pj are algebraic numbers. Thus the problem of 
reaching a hyperplane infinitely often reduces to the Infinite Zeros 
Problem', given an exponential polynomial /, decide whether / has 
infinitely many non-negative real zeros. Note that since / is an an¬ 
alytic function on the whole real line it only has finitely many zeros 
in a bounded interval. Thus the Infinite Zeros Problem is equivalent 
to asking whether the set of zeros of / is unbounded. 

This paper is concerned with the decidability of the Infinite 
Zeros Problem. In order to formulate our main results, recall that 
exponential polynomials can equivalently be characterised as the 
solutions of ordinary differential equations 

-b Un-i/*'" -b . . . -b aof = 0 , (1) 

with the coefficients Oj and the initial conditions (0) being real 
algebraic numbers for j G {0,..., n — 1}. We say that / has order 
n if it satisfies a linear differential equation of the form (TJ. 

Our main results concern both decision procedures and hard¬ 
ness results for the Infinite Zeros Problem. We show that the prob¬ 
lem is decidable for exponential polynomials of order at most 7. 
With regards to hardness, we exhibit a reduction to show that de¬ 
cidability of the Continuous Infinite Zeros Problem for instances 
of order at least 9 would entail major advancements in the field of 
Diophantine Approximation, namely the computability of the La¬ 
grange constants of arbitrary real algebraic numbers. 

Let us expand on the significance of the above hardness result. 
Essentially nothing is known about the Lagrange constant of any 
real algebraic number of degree three or above. For example, it has 
been a longstanding open problem since the 1930s whether some 
real algebraic number of degree at least three has strictly positive 
Lagrange constant and, on the other hand , whether s ome such num¬ 
ber has Lagrange constant 0 (see, e.g., ( lGuvll20()^ ). These ques¬ 
tions are often formulated in terms of the simple continued fraction 
expansion of a real number a, which has unbounded elements if 
and only if a has Lagrange constant 0. 

The reader will notice that there is a gap between our decid¬ 
ability and hardness results for exponential polynomials of order 8. 
We claim decidability in this case but defer the details to a longer 
version of this paper. 

Another way to calibrate the difficulty of the Infinite Zeros 
Problem for an exponential polynomial f{t) = 
is in terms of the dimension of the Q-vector space spanned by 
{Im(Q:j) : j — 1,..., k}. We show decidability in case this space 
is one-dimensional and we observe that the above hardness result 
already applies in the two-dimensional case. 
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1.1 Related Work 

Closely related to the Infinite Zeros Problem is the problem of 
whether an exponen tial polynomial has some zero. This prob¬ 
lem is considered in dBell et al.ll2010h under the name Continuous 
Skolem-Pisot Problem. The techniques considered in the present 
paper are relevant to the latter problem, but significant extra diffi¬ 
culties arise in this new setting since we cannot discount the be¬ 
haviour of / on some bounded initial segment of the reals. Our 
work on the Continuous Skolem-Pisot Problem will be reported 
elsewhere. 

There is a natural discrete analog of the Infinite Zeros Problem: 
given a linear recurrence sequence, determine whether it has in¬ 
finitely many zero terms. The decidab ility of the latter problem wa s 
established by Berstel and Mignotte dBerstel and Mignottelll97dl . 
The problem of deciding whether a given linear recurrence se¬ 
quence has some zero term is called Skolem’s Problem. This is a 
longstanding and celebrated open problem which essentially asks 
to give an effective proof of the Skolem-Mahler-Le ch Theore m for 
linear recurrences; see, e.g., the exposition of Tao dTadl200^ Sec¬ 
tion 3.9). 

Macintyre and Wilkie dMacintvre and Wilkie! 199^1 showed de¬ 
cidability of the first-order theory of (R,-f, x,0, l,<,exp), the 
real field with exponentiation, subject to Schanuel’s Conjecture in 
transcendence theory. In this paper we are concerned with the com¬ 
plex exponential function, and we do not use this result. Moreover, 
although we do make use of transcendence theory, all the results in 
this paper are unconditional. 

2. Mathematical Background 

2.1 General Form of a Solution 

We recall some facts about the general form of solutions of ordi¬ 
nary linear differential equations. Consider a homogeneous linear 
differential equation 

/(") + + ... + Co/ = 0 (2) 

of order n. The characteristic polynomial of (Ell is 

xix) := -I- -I- ... -I- Co . 

If A is a root of x of multiplicity m, then the function f{t) = 
satisfies El for / = 0,1,..., m — 1. There are n distinct linearly 
independent solutions of El having this form, and these span the 
space of all solutions. 

Let the distinct roots of x be Ai,..., Afc, with respective multi¬ 
plicities mi,..., mfc. We refer to Ai,..., Aj, as the characteristic 
roots of the differential equation. We also refer to the characteris¬ 
tic roots of maximum real part as dominant. Write Xj = Vj -f ioj 
for real algebraic numbers rj,aj, j = 1,..., fc. It follows from 
the discussion above that, given real algebraic initial values of 
f(0), f'{0 ),..., the uniquely defined solution / of (El 

can be written in one of the following three equivalent forms. 

1. As an exponential polynomial 

k 

3 = 1 

where each Pj is a polynomial with (complex) algebraic coef¬ 
ficients and degree at most rrij — 1. 

2. As a function of the form 

k 

/(f) = ^e^\Pj{t)cos{ajt) + Qj{t)s\a{ajt)) 
i=i 

where the polynomials Pj , Qj have real algebraic coefficients 

and degrees at most rrij — 1. 


3. As a function of the form 

k mi —1 

/(f) = E E cos{ajt + pjj) 

1 = 1 i=0 

where bjj is real algebraic and algebraic for each j, 1. 

We refer the reader to dBell et al.ll2010l Theorem 7) for details. 

2.2 Number-theoretic tools 

Throughout this paper we denote by A the set of algebraic numbers. 
Recall that a standard way to represent an algebraic number a 
is by its minimal polynomial M and a numerical approximation 
of sufficient acc uracy to distinguish a from the other roots of 
M dCohenI[19931 Section 4.2.1). Given two algebraic numbers a 
and /3 under this representation, the Field Membership Problem is 
to determine whether jd € Qja) and, if so, to return a polynomial 
P with rational coefficients such that jd — P(a). This problem 
can be decided using the LLL algorithm, see dCohen|[T993l . Section 
4.5.4). 

Given the characteristic polynomial x of a. linear differential 
equation we can compute approximations to each of its root s 
Ai,..., An to within an arbitrarily small additive error dPanll996h . 
Moreover, by repeatedly using an algorithm for the Field Member¬ 
ship Problem we can compute a primitive element 0 for the splitting 
field of X and representations of Ai,..., An as polynomials in 8. 
Thereby we can determine maximal (Q-linearly independent sub¬ 
sets of {Re(Aj) '■ 1 < j < n} and {Im(Aj) : 1 < / < n}. 

We now move to some techniques from Transcendental Num¬ 
ber Theory on which our results depend in a critical way. The 
following theo r em was originally proven in 1 934 by A. Gelfond 
dGelfondi 19341 ; iGelfond and yinogradov|[T934h and independently 
by T. Schneider dSchneidej Il935al lbll. settling Hilbert’s seventh 
problem in the affirmative. 

Theorem 1. ( Gelfond-Schneider) If a and b are algebraic numbers 
with a yt 0,1 and b ^ Q, then a*’ is transcendental. 

The following lemma, proven in dBell et al.j |2()jT ), is a useful 
consequence of the powerful Baker’s Theorem dBakei 119751 Theo¬ 
rem 3.1): 

Lemma 2. iBell etal.\\201(X Lemma 13) Let a,b G R n A foe 
linearly independent over ((J and let tpi, ip 2 be logarithms of al¬ 
gebraic numbers, that is, € A. There exist effective 

constants C,N,T > 0 such that for all t > T, at least one of 
1 — cos{at d- ipi) > Cjt^ and 1 — cos{bt -|- 952 ) > Cholds. 

Another necessary tool is a version of Kronecker’s well-known 
Theorem in Diophantine Approximation. 

Theorem 3. (Kronecker, appears in ItHardy and WriehAl99f\) ) Let 
Ai,..., Am and xi,..., Xm be real numbers. Suppose that for all 
integers u\,..., Um such that uiXi + • • • -I- UmXm G Z, we also 
have uixi -1- • • • -I- UmXm G Z, that is, all integer relations among 
the Xj also hold among the Xj (modulo Z). Then for all e > 0, 
there exist p G Z™' and n G N such that \nXj — Xj — Pj\ < e 
for all 1 < j < m. In particular, if 1, Xi,. . . , Am are linearly 
independent over Z, then there exist such n G N and p G Z"* for 
all X G R'" and e > 0. 

A direct consequence is the following: 

Lemma 4. Let ai,... , Om G R n A foe linearly independent 
over ((J and let ip\,..., prn. G R. Write x mod 27r to denote 
minfcgz \x -I- 2A:7r| for any x G R. Then the image of the mapping 
h(f) : R>o —>■ [0, 27r)”^ given by 

h{f) = {{a\t -I- tpi) mod 2it,..., {omt <pm) mod 27r) 


2 


2016/5/10 





























is dense in [0, 27r)"‘. Moreover, the set 

{h{t) I (ait + ipv) mod 27r = 0} 
is dense in {0} x [0, 27r)™“^. 

Proof. For the first part of the claim, note that the linear inde¬ 
pendence of l,ai/27r,... ,am/27r follows from the linear inde¬ 
pendence of ai,...,Om and the transcendence of tt. Then hy 
Kronecker’s Theorem, the restriction {h{t) 11 G N} is dense in 
[0, 271)"*, so certainly the whole image of h{t) must also he dense 
in [0, 27r)™'. For the second part, the trajectory h{t) has zero first 
coordinate precisely when t = —ipijai -I- 2mt for some n G Z, at 
which times the trajectory is 

g{n) h ^^ -f 2n7r^ 

= {0} X ^ 

V “1 

As before, we have that {1, , 27ram,/ai} are linearly 

independent over Q from the linear independence of ai,, Om 
and the transcendence of tt, so applying Kronecker’s Theorem to 
the last m — 1 components of this discrete trajectory yields the 
second part of the claim. □ 

2.3 First-Order Theory of the Reals 

We denote by C the first-order language (R, -I-, x, 0, 1, <). Atomic 
formulas in this language are of the form P{xi,..., Xn) = 0 and 
P{xi,..., Xn) > 0 for P G Z[xi,..., Xn] a polynomial with 
integer coefficients. A set X C R" is definable in C. if there exists 
some P-formula 4>{x) with free variables x which holds precisely 
for valuations in X. Analogously, a function is definable if its graph 
is a definable set. 

We denote by T'/i(R) the first-order theory of the reals, that 
is, the set of all valid sentences in the language C. It is worth re¬ 
marking that any real algebraic number is readily definable within 
C using its minimal polynomial and a rational approximation to 
distinguish it from the other roots. Thus, we can treat real alge¬ 
braic numbers constants as built into the language and use them 
freely i n the constru ction of formulas. A celebrated result due to 
Tarski ( lTarskj|ll95 Ih is that the first-order theory of the reals ad¬ 
mits quantifier elimination: that each formula (x) in C is equiv¬ 
alent to some effectively computable formula (f> 2 {x) which uses 
no quantifiers. This immediately entails the decidability of Th{M.). 
It also follows that sets definable in £. are precisely the semialge- 
braic sets. Tarski’s original result had non-elementary complexity, 
but impr ovements follow ed, culminating in the detailed analysis of 
Renegar iRene g^l 199211 . 

Decidability and geometrical properties of definable sets in the 
first-order theory of the structure Cexp = (R, -f, x, 0,1, <, exp), 
the reals with exponentiation, have been exp lored by a number of 
authors. Most notably, Wilkie (IWilkielll99al showed that the the¬ 
ory is o-minimal and Macintyre and Wilkie dMacintvre and Wilkiel 
Il996ll showed that if Schanuel’s conjecture is true then the theory 
is decidable. We will not need the above two results in this paper, 
however we use the following, which is very straightforward to es¬ 
tablish directly. 

Proposition 5. There is a procedure that, given a semi-algebraic 
set S C R*’ and real algebraic numbers ai,..., ak, returns an 
integer T such that {f > 0 : ,. ■., 6“*"*) G 5} either contains 

the interval {T, oo) or is disjoint from (T, oo). The procedure also 
decides which of these two eventualities is the case. 

Proof. Consider a polynomial P G Z[ui,..., ttfc]. For suit¬ 
ably large t the sign of ..., is identical to the 


sign of the coefficient of the dominant term in the expansion of 
P(e“'^*,..., 6“*^*) as an exponential polynomial. It follows that 
the sign of P(e“^*,.. ., 6“*=*) is eventually constant. It is more¬ 
over clear that one can effectively compute a threshold beyond 
which the sign of P(e“^*,..., e“'=*) remains the same. Since 
the set S is defined by a Boolean combination of inequalities 
P(mi, .. ., Uk) ~ 0, for ~ G {<, =}, the proposition immedi¬ 
ately follows. □ 

2.4 Useful Results About Exponential Polynomials 

We restate two useful theorems due to Bell et al. dBell et akllTOlOh . 

Theorem 6. iBell et al.l201(X Theorem 12} Exponential polynomi¬ 
als with no real dominant characteristic roots have infinitely many 
zeros. 

Theorem 7. iBell et a/1l20iCl Theorem 15) Suppose we are given 
an exponential polynomial whose dominant characteristic roots are 
simple, at least four in number and have imaginary parts linearly 
independent over Q. Then the existence of infinitely many zeros is 
decidable. 

3. One Linearly Independent Oscillation 

In this section we consider exponential polynomials f{t) — 
Pj{t)e^^^ under the assumption that the span of {Im(Aj) : 
y = 1,..., A:} is a one-dimensional Q-vector space. In this case we 
can use fundamental geometric properties of semi-algebraic sets to 
decide whether or not / has finitely many zeros. 

Theorem 8. Let f{t) = Pj(t)e^P be an exponential poly¬ 

nomial such that the span of {Im(Aj) \ j = 1,... ,lt} is a one¬ 
dimensional ij-vector space. Then the existence of infinitely many 
zeros of f is decidable. 

Proof. Write \j = Oj -|- ibj, where aj, bj are real algebraic num¬ 
bers for j = 1,..., fc. By assumption there is a single real algebraic 
number b such that each bj is an integer multiple of b. Recall that 
for each integer n, both cos{nbt) and sm{nbt) can be written as 
polynomials in sm{bt) and cos{bt) with integer coefficients. Using 
this fact we can write / in the form 

fit) = Qit, ..., cos(6f), sin(6 f)), 

for some multivariate polynomial Q with algebraic coefficients. 
Now consider the semi-algebraic set 

P:={(n,s)GR''+^Q(uo,...,w,,^,^) = o} . 

Recall that | ^ j : s G r| comprises all points in the 

unit circle in R^ except (—1,0). Indeed, given 9 G {—it,it), 
setting s := tan(0/2) we have cos(9) — and sin(0) = 

It follows that f(t) = 0 and cos(6f) yt —1 imply that 
(f,e“i*,...,e“*^*,tan(bf/2)) G E. 

By the Ce ll Decomposition Theorem for semi-algebraic sets 
iMarkeil2002h . there are semi-algebraic sets Ci,... , Cm C 
Pi,..., Dm C and continuous semi-algebraic functions 

: Pj —t R such that E can be written as a disjoint 
union P = Ci U ... U Cm, where either 

Cj = {(li, s) G : u G Dj A s = Cj('w)} (3) 

or 

Cj = {{u, s) G R''+" -.uGDjA (n) < s < cf ^ (n)} (4) 


<piaj 


mod 27r 


2<j<m 
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Figure 1. Intersection points of r]j {t) and tan(6f/2). 


Moreover such a decomposition is computable from E. Clearly 
then 

m 

{f £ R : f{t) = 0} C y G R : (t, ..., G Dj}\jZ , 

where Z := {f £ R : cos(6f) = —1}. 

The restriction of the exponential polynomial / to Z is given 
by /(f) = Q(t, ..., —1,0). Since this expression is a 

linear combination of terms of the form for real algebraic r, 
for sufficiently large t the sign of f{t) is determined by the sign 
of the coefficient of the dominant term. Thus / is either identically 
zero on Z (in which case / has infinitely many zeros) or there exists 
some threshold T such that all zeros of / in Z lie in the interval 
[0,T]. 

We now consider zeros of / that do not lie in Z. There are two 
cases. First suppose that each set {f £ R : {t, ..., e“''*) £ 

Dj} is bounded for j = 1,..., m. In this situation, using Propo¬ 
sition |5] we can obtain an upper bound T such that if f{t) = 
0 then t < T. On the other hand, if some set {t £ R>o : 

... ,e“''*) £ Dj} is unbounded then, by Proposition!^ 
it contains an infinite interval (T, oo). We claim that in this case / 
must have infinitely many zeros t > 0. We first give the argument 
in the case Cj satisfies (gil. 

Define rij{t) = ..., 6“*=*) for t £ (T, oo). Then for 

t £ (T, oo) \ Z, 

/(f) = 0 <= (f, ..., e“*'*,tan(6f/2)) £ Cj 

(t, , 6°’’°^) £ Dj A rij(t) = tan(6f/2). 

In other words, / has a zero at each point t £ (T, oo) \ Z at 
which the graph of rjj intersects the graph of tan(6f/2). Since 
rjj is continuous there are clearly infinitely many such intersection 
points, see Figure[T] 

The case when Cj satisifes © is handled similarly. In fact, 
this case cannot arise at all, since by the above argument, if Cj 
satisfies then / has a non-trivial interval of zeros. This is 
impossible, since / is analytic, and hence has only isolated zeros. 
This completes the proof. □ 


4. Decidability Results up to Order 7 

We now shift our attention to instances of the Infinite Zeros Prob¬ 
lem of low order. In particular, for exponential polynomials corre¬ 
sponding to differential equations of order at most 7, we establish 
decidability of the Infinite Zeros Problem. 

Theorem 9. The Infinite Zeros Problem is decidable for differential 
equations of order at most 7. 

Proof. Suppose we are given an exponential polynomial / of order 
at most 7. Sort the characteristic roots according to their real parts, 
and let rj denote throughout the y-th largest real part of a character¬ 
istic root. We will refer to the characteristic roots of maximum real 


part as the dominant characteristic roots. Let also mul(\) denote 
the multiplicity of A as a root of the characteristic polynomial of /. 

We will now perform a case analysis on the number of dominant 
characteristic roots. By Theorem!^ it is sufficient to confine our 
attention to exponential polynomials with an odd number of dom¬ 
inant characteristic roots. Throughout, we rely on known general 
forms of solutions to ordinary linear differential equations, outlined 
in SectionjzTj 

Case I. Suppose first that there is only one dominant, necessarily 
real, root r. Then if we divide / by we have: 

^=Pi(f)+o(e(’'-’'>‘), 

as the contribution of the non-dominant roots shrinks exponentially, 
relative to that of the dominant root. Thus, for large f > 0, the sign 
of fit) matches the sign of the leading coefficient of Pi{t), so / 
cannot have infinitely many zeros. 

Case II. We now move to the case of three dominant character¬ 
istic roots: r and r ± ia, so that 

= Pi(f) -b P 2 (f) cos(of) + P3{t) sin(of) -b O , 

where Pi, P 2 , Ps G (R fl A)[a;] have degrees di deg(Pi) < 

def 

mul{r) — 1 and da = deg(P 2 ) = deg(P 3 ) < mul(r ± ai). 

Case Ila. Suppose di > da. Now, it is easy to see that for large 
t the sign of f{t) matches the sign of the leading coefficient pi of 

Pi: 

so clearly some bound T exists such that t > T /(f) 7 ^ 0. 
Similarly, if da > di, then /(f) clearly has infinitely many zeros. 
Indeed, if pa, P 3 are the leading coefficients of Pa, P 3 , respectively, 
then we have: 

=P 2 Cos(af) -bp 3 sin(af) -bC>(l/f) +0 

= + V" ) 0 ( 1 /^) + 0 ( 


where ip £ [0, 27r) with tan((b>) = —pa/pa, so / is infinitely often 
positive and infinitely often negative. 

Thus, we can now assume di = da. Notice that since the order 
of our exponential polynomial is no greater than 7, we must have 
di = da ^ 2. 

Case Ilb. Suppose that di = da = 2. Then our function is of 
the form 

= f(blcos(af-bv5i)+P)+(C'cos(af-bv52)+P)-be*'''^~’’^*P, 
e’’* 

for constants A, B, C, D, F,a G RnA with a > 0 and 6*“^^, £ 

A. In this case, Theorem|9|follows from LemmafT^in Section iTAl 

Case lie. Suppose that di = da = 1, so that 

= Ai cos(af -b b^i) -b Aa -b e^’’^“’''*Pi(f), 

where Ai,A 2 ,a £ R n A, a > 0, 6*“^^ £ A and Pi is an 
exponential polynomial with dominant characteristic root whose 
real part is 0. Consider first the magnitudes of Ai and A^. If 
I All > jAaj, then the term Ai cos(af -b v^i) makes / change 
sign infinitely often, so / must have infinitely many zeros. On the 
other hand, if |Ai| < jAaj, then / is clearly ultimately positive 
or ultimately negative, depending on the sign of Aa. The remaining 
case is that | Ai| = lAaj. Dividing / by Aa, replacing (pi by pi -bTr 
if needed and scaling constants by Aa as necessary, we can assume 
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the function has the form; 


/w 

grt 


1 — cos(af + (^i) + 


We now enumerate the possibilities for the dominant characteristic 
roots of the exponential polynomial Fi, that is, the characteristic 
roots of / with second-largest real part. Since / has order at most 
7, there are the following cases to consider: 


• Fi has four simple, necessarily complex, dominant roots, so 
that 


M = i 

ert 


cos(af + ifii) 


_l_ g (’'2 ’’)*(5cos(6f -I- (^ 2 ) -I- Ccos(cf -I- ifa)), 


where A,B,C,a,b G K n A, a, 6 > 0, 6*“^^, G A and F 
is an exponential polynomial of order at most 2 whose dominant 
characteristic roots have real part equal to 0. In this case. Theorem 
[^follows from Lemma [T4l in Section |4^ 

Case IV. Finally, suppose there are seven dominant characteris¬ 
tic roots: r, r±ai, r±bi and r±ci. Since we are limiting ourselves 
to instances of order 7, these roots must all be simple, and there can 
be no other characteristic roots. Thus, the exponential polynomial 
has the form 

= Acosiat + (fi) + B cosfbf -I- 1 P 2 ) + C cosict + ips) + D, 

grt 

with A, B,C, D, a,b,c G R n A, a, b, c > 0 and , ■ • •, £ 

A. In this case, Theoreml^follows from Lemma [TS] in Section l43] 

□ 


where B, C, b, c G R IT A with b, c > 0 and G A. In 

this case, Theorem|^follows from LemmafT^in Section ldTI 
• Fi has some subset of one real and two complex numbers as 
dominant roots, all simple, so that 


= 1 - cos(af -I- ifi) 

_l_ g(’’2-’')‘(5cos(bf -F '^ 2 ) -F C) -F e^’’^“’''‘F 2 (f), 


where B,C,b G R T A, b > 0, G A and F 2 is an ex¬ 
ponential polynomial with dominant characteristic root whose 
real part is 0. In this case. Theoremfollows from Lemma fTTI 
in Section IdAl 

• Fi has a repeated real and possibly two simple complex domi¬ 
nant roots, so that 

= 1 - cos(af -F ipi) 
e 

+ e^’'^“’'^‘(Bcos(bf -F ¥^ 2 ) -F P{t)) + 


where B,bGRnA, b>0, G A, and P G (R T A) [x] is 
non-constant. Now, if the leading coefficient of P is negative, 
then / will be infinitely often negative (consider large times 
t such that cos{at -F v^i) = 1) and infinitely often positive 
(consider large times t such that cos{at -\- tpi) = 0 ), so / must 
have infinitely many zeros. On the other hand, if the leading 
coefficient of P is positive, then it is easy to see that / is 
ultimately positive. 

• Pi has a repeated pair of complex roots, so that 
^ = l-cos(af-F¥3i) 

grt 

_l_ g,P^~A^{^Btcos{bt -F (p 2 ) + C cos(bf -F T’s)), 


where PjC, bGROA, b>0 and G A. In this case, 

Theoreml^follows from LemmafT^in Section l40l 


Case III. We now consider the case of five dominant charac¬ 
teristic roots. Let these be r, r ± ai and r ± bi. If r ± ai are 
repeated, i.e., mul{r ± ai) > 2 , then we must have mul{r) = 
mul{r ± bi) = 1 , since otherwise the order of our exponential 
polynomial exceeds 7. Then by an argument analogous to Case Ila 
above, / must have infinitely many zeros. The situation is sym¬ 
metric when mul{r ± bi) > 2. Similarly, if mul{r) > 2, then 
mul{r ± ai) = mul{r ± hi) = 1 , since otherwise the instance 
exceeds order 7. Then by the same argument as in Case Ila, f is ul¬ 
timately positive or ultimately negative. Thus, we may assume that 
all the dominant roots are simple, so the exponential polynomial is 
of the form: 

= Acos{at + 1 ^ 1 ) + Bcos{bt + ^ 2 ) -F C -F e’'’'^“'^^*P(f), 

grt 


In the remainder of this section, we provide the technical lem¬ 
mas invoked throughout the proof of Theorem[3 

4.1 One dominant oscillation 

Lemma 10. Let A, B,a,b,r G R T A where a,b,r > 0. Let 
, 7^2 £ R be such that G A. Suppose also that a, b are 

linearly dependent over and that whenever l — coa{at+tpi) = 0, 
it holds that A cos(bf -F 732 ) + P > 0. Define the function 

f{t) = 1 — cos{at + (fi) + e~^*{Acos{bt + ip 2 ) -F B). 

Then f{t) = that is, there exist effective constants T > 0 

and c > 0 such that for t > T, we have f{t) > 


Proof. The case of A = 0 is easy: by the premise of the Lemma, we 
have B > 0 and then f{t) > Be~^* for all t. Thus, assume A 7 ^ 0 
throughout. Let the linear dependence between a, b be given by 
ani — bn 2 — 0 for ni,n 2 G N coprime and letC be the equivalence 
class of —pija modulo 2n/a, that is. 


Q def J —7^1 + 2fc7r 


k G 


We will refer to C as the set of critical points throughout. 

It is clear that at critical points, we have 1 — cos(af -F T^i) = 0. 
Moreover, the linear dependence of a, b entails that for each fixed 
value of (cos(af), sin(of)), there are only finitely many possible 
values for (cos(bf), sin(bf)). Indeed, we have 

e”*"* G I LO an n 2 -th root of unity}, 


so in particular, for t G C,we have 

e**** G I uj an n 2 -th root of unity}. 

Thus, the possible values of (cos(bf), sin(bf)) for t critical are 

def 

algebraic and effectively computable. Let M — min{Acos(bf + 
7 > 2 ) + P|f£C}.By the premise of the Lemma, we have M > 0. 

Let ti,t 2 ,... ,tj,. .. be the non-negative critical points. Note 
that by construction we have \tj — tj-i\ = 2 tx fa. For each tj, 
define the critical region to be the interval [tj — 5, tj + b], where 

- def M 
~ 2|A|b' 


Let g{t) = Acos(bf -F P 2 ) + B and notice that g'{f) < \A\h 
everywhere. We first prove the claim for t inside critical regions: 
suppose t lies in a critical region and let j minimise \t — tj\ < 5. 
Then by the Mean Value Theorem, we have 

\9it)-g{ij)\ < |f-fj||A|b<<5|A|b= y, 


/X , ^ M M 
9{t) > fffe) “ Y - Y’ 


5 


2016/5/10 





whence f{t) > e-^^g{t) > Me-''V 2 = 

Now suppose t is outside all critical regions and let j minimise 
\t — tj\. Since the distance between critical points is 2'K/a by 
construction, we have a\t — tj \ < tt. Therefore, 

1 — cos{at + ifii) = 1 — cos{at — atj) > 2 ^^^^ 

(aS)^ a^M'^ 

> -—— = - > 0 

2 8 |A| 2 fe 2 

Thus, there exists a computable constant D > 0 such that f{t) = 
1 — cos{at + 9 P 1 ) + e~'~^g{t) > D for all large enough t outside 
critical regions. 

Combining the two results, we have f{t) = every¬ 
where. □ 

Lemma 11. Let C, D, a, b, ri, r 2 be real algebraic numbers such 
that a,h,r\,r 2 > 0 and C, D are not both 0. Let also </ 5 i, V 32 € M 
be such that € A. Define the exponential polynomial f 

by 

f{t) = 1 — cos(af -I- pi) 

+ e""i*(Ccos(6f + p2) + D)+ 

Here F is an exponential polynomial whose dominant character¬ 
istic roots are purely imaginary. Suppose also that / has order at 
most 7. Then it is decidable whether f has infinitely many zeros. 

Proof. Notice that the dominant term of / is always non-negative, 
so the function is positive for arbitrarily large t. Thus, /(f) = 0 for 
some t if and only if f(t) < 0 for some t, and analogously, / has 
infinitely many zeros if and only if f{t) < 0 infinitely often. We 
can eliminate the case \D\ > \C\, since then / is clearly ultimately 
positive or oscillating, depending on the sign of D. Thus, we can 
assume |I 1 | < IC]. 

We now consider two cases, depending on whether a/b £ Q. 
Case I. Suppose first that a, b are linearly independent over Q. 
By Lemma|4l the trajectory {at -\- pi mod 27r, bt p 2 mod 27r) 
is dense in [ 0 , 27r)^, and moreover the restriction of this trajectory 
to at -b pi mod 27r = 0 is dense in {0} x [0, 27r). 

If |D| < I Cl, then we argue that / is infinitely often negative, 
and hence has infinitely many zeros. Indeed, |D| < |C| entails the 
existence of a non-trivial interval I C [0, 27r) such that 

t mod 27r £ 7 C cos{bt -b P 2 ) -b 7) < 0. 

What is more, we can in fact find e > 0 and a subinterval I' C I 
such that 

t mod 27r £ 7 ' C cos{bt + P 2 ) + D < —e. 

Thus, by density, 1 — cos(af + pi) =0 and C cos(6f -b P 2 ) -b 77 < 
—e will infinitely often hold simultaneously. Then just take t large 
enough to ensure, say, |e“’'^*7^(f)| < e/2 at these infinitely many 
points, and the claim follows. 

Thus, suppose now |C| = |77|. Replacing p 2 hy p 2 -\- tt if 
necessary, we can write the function as; 

/(f) = 1 — cos(af -b pi) 

+ 77e""i‘(l - cos(6f -b P 2 )) + e-<’'^+’'=^‘T’(f). 

As a,b are linearly independent, for all f large enough, 1 — 
cos(af + pi) and 1 — cos(6f -b P 2 ) cannot simultaneously be ‘too 
small’. More precisely, by Lemmaj^ there exist effective constants 
E,T,N >0 such that for all t >T, we have 

1 — cos(af -b pi) > E/t^ or 1 — cos(6f -b P 2 ) > E/t^. 


Now, if 77 < 0, it is easy to show that / has infinitely many zeros. 
Indeed, consider the times f where the dominant term 1 — cos(af -b 
pi) vanishes. For all large enough such f, since t~^ shrinks more 
slowly than we will have 

/(f) = e~’'’^*77(l — cos(bf -b P 2 )) + F{f) 

< e~''^\EDt-^ + e“’'"‘T’(f)) 

< e~^^^-EDt-^ 

2 

<0, 

so / has infinitely many zeros. Similarly, if 77 > 0, we can show 
that / is ultimately positive. Indeed, for all f large enough, we have 

/(f) > e~’'’^*77(l — cos(bf -b <^> 2 )) -b F{f) 

> e-'^^^DEt-^ + e-(''i+"=^‘7’(f) 

>0, 

or 

/(f) > 1 - cos(af -b pi) + e"^"i+"")‘T’(f) 

> Et~^ -b e"^"i+"">*T’(f) 

> 0. 

Therefore, / has only finitely many zeros. 

Case II. Now suppose a, b are linearly dependent. By the 
premise of the Lemma, the order of F is at most 2 (in fact, at 
most 1 if 77 yf 0). However, by Theorem]^ the claim follows im¬ 
mediately for all cases in which the characteristic roots of F are 
all real or complex but with frequencies linearly dependent on a. 
Thus, the only remaining case to consider is the function 

/(f) = 1 — cos(af -b <b>i) 

-b cos{bt -b P 2 ) + H cos(cf -b pf), 

where 77, c£]RnA, c>0 and a/c ^ Q. 

As explained at the beginning of the proof of LemmafTOl due to 
the linear dependence of a, b over Q, when 1 — cos(af -b ‘b’l) = 0, 
there are only finitely many possibilities for the value of C cos(6f-b 
P 2 ), each algebraic, effectively computable and occurring period¬ 
ically. If at least one of these values is non-positive, then by the 
linear independence of a, c over Q, we will simultaneously have 
1 —cos(af-b</5i) — 0, Ccos{bt+p 2 ) < 0 and 77 cos(cf-bv33) <0 
infinitely often, which yields /(f) < 0 infinitely often and en¬ 
tails the existence of infinitely many zeros. On the other hand, 
if at the critical points 1 — cos(af + pi) = 0 we always have 
C cos(6f -\- P 2 ) > Q, then by Lemma [Tol we have 

1 — cos(of -b pi) -b cos(6f -b P 2 ) = 0(e~'’’-*), 

whereas obviously 

|e-("i+’'=)*77cos(cf -b ¥^3)| = C7(e-("i+"^>*). 

If follows that / is ultimately positive and hence has only finitely 
many zeros. □ 

Lemma 12. Let A, B, a,b,c,r be real algebraic numbers such 
that a, b,c,r > 0, A, B yf 0. Let also yji, 1 /P 2 , (^3 C M. be such 
that 6*“^^, £ A. Define the exponential polynomial f by 

f{t) = l — cos{ct-\-p 3 )-\-e~^^{Acos{at-\-pi)-\-Bcos{bt-\-p 2 )). 

Then it is decidable whether / has infinitely many zeros. 
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Proof. We argue the function is infinitely often positive and in¬ 
finitely often negative by looking at the values of t for which the 
dominant term 1 — cos(cf + ips) vanishes. This happens precisely 
at the times t = —{tpzP 2kn )/c for A: € Z, giving rise to a discrete 
restriction of /: 


9{k) 


B cos I k 


2'iTb bipz 


+ pi 


This is a linear recurrence sequence over R of order 4, with char¬ 
acteristic roots and In particular, it has no 

real dominant characteristic root. It is well-known that real-valued 
linear recurrence sequences with no dominant real characteristic 
root are infi nitely often positive an d infinitely often negative: see 
for example dGvori and Ladasri99lL Theorem 7.1.1). Therefore, by 
continuity, / must have infinitely many zeros. □ 


Lemma 13. Let A, B, a, b, r be real algebraic numbers such that 
a,b,r > 0, A ^ 0. Let also (fii,p 2 ,P 3 G R foe such that 
g»¥>i ^ g*<e 2 ^ g »¥’3 ^ Define the exponential polynomial f by 

f{t) = l — cos{at+(pi) + e~''*{Atcos{bt+(p 2 ) + Bcos{bt+p 3 )). 

Then it is decidable whether f has infinitely many zeros. 


Proof. If a/fo G Q, then the claim follows immediately from 
Theorem[ 8 l If a/fo ^ Q, then by Lemma|4l it will happen infinitely 
often that 1 — cos(af -|- pfi) — 0 and Atcos{bt -|- P 2 ) < —\A\t/2. 
Then clearly f(t) < 0 infinitely often. Since /(f) > 0 infinitely 
often as well, due to the non-negative dominant term 1 — cos (at + 
1 ^ 1 ), it follows that / has infinitely many zeros. □ 

4.2 Two dominant oscillations 

Lemma 14. Let A, B, C, a, b, r be real algebraic numbers such 
that a,b,r > 0, a b and A,B,C 7 ^ 0. Let also pi,p 2 G R foe 
such that G A. Define the exponential polynomial f by 

f{t) = Acos{at + pi) + B cos(fof -I- P 2 ) -I- C -I- e~'"*F{t). 

where F is an exponential polynomial whose dominant character¬ 
istic roots are purely imaginary. Suppose also f has order at most 
8 . It is decidable whether f has infinitely many zeros. 


We are now interested in the extrema of P(sin(f), cos(f)). Let 

def 

Ml = min P{x,y) = 

X‘^+y^ = l t>U 

def 

M 2 = max P(x,v) — maxa(f). 

i>o 

We can construct defining formulas <j>i{u), 1 / 2 (u) in the first-order 
language €. of real closed fields for Mi, M 2 , so that each 
holds precisely for the valuation u = Mj . Then performing quan- 
tifier eliminati on on these formulas using Renegar’s algorithm 
dRenegaj|l992h . we convert cj>i, (f >2 into the form 

(j}j{u) = \f f\Pi^k{u) 0 , 

i k 

where Pi^k are polynomials with integer coefficients and each 
is either < or =. Now ifij (u) must have a satisfiable disjunct. Using 
the decidability of the theory rfo(R), we can readily identify this 
disjunct. Moreover, since (u) has a unique satisfying valuation, 
namely u — Mj, this disjunct must contain at least one equal¬ 
ity predicate. It follows immediately that Mi, M 2 are algebraic. 
Moreover, we can effectively compute from fijiu) a representa¬ 
tion for Mj consisting of its minimal polynomial and a sufficiently 
accurate rational approximation to distinguish Mj from its Galois 
conjugates. By an analogous argument, the pairs (sin(f),cos(f)) 
at which P(sin(f), cos(f)) achieves the extrema Mi, M 2 are also 
algebraic and effectively computable. 

We now perform a case analysis on the signs of Mi and M 2 . 

• First, if 0 < Mi < M 2 , then f(t) cannot have infinitely many 
zeros: if t is large enough to ensure |/3(f)| < Mi, we have 
/(f) > 0. 

• Second, if Mi < M 2 < 0, then by the same reasoning, the 
function will ultimately be strictly negative. 

• Third, if Mi < 0 < M 2 , then / oscillates around 0: for 
all t such that a(i) = Mi < 0 and large enough to ensure 
|/?(f)| < |Mi|, we will have /(f) < 0, and similarly, for large 
enough f such that a(f) = M 2 > 0 , we will have /(f) > 0 , so 
the function must have infinitely many zeros. 

• Next, we argue that the case Mi = M 2 = 0 is impossible. 
Indeed, if Mi = M 2 = 0, then Q:(f) = P(sin(f), cos(f)) is 
identically zero, and the same holds for all derivatives of a(f). 
Thus, from 01 {t) = of" (t) = 0, we have 

0 = —Am sin(mf -\- pi) — Bn sin(nf -F P 2 ), 

0 = Am^ sin(mf -F pi) -F Bn?' sin(nf -F P 2 ). 


Proof. If the frequencies a, fo of the dominant term’s oscillations are 
linearly independent over Q, then the claim follows immediately by 
Theorem|7] Therefore, assume na — mb = 0 for some n,m £ N'*'. 
Notice that a b guarantees n m. We perform the change of 
variable f —^ fm/a, so that: 

/(f) = A cos(mf-FvJi)-FB cos(nf-Fv52)+C'-Fe~^"‘*^“P(fm/a). 

Using the standard trigonometric identities, we express the domi¬ 
nant term as a polynomial in sin(f), cos(f): 

/(f) = P(sin(f), cos(f)) -\- e~™*'^°'F{tm/a), 

where P G (R IT A)[x, y] has effectively computable coefficients. 
It is clear that the dominant term is periodic. It is immediate from 
the definition of exponential polynomials and the premise of the 
Lemma that F{tm/a) = Fzft) is an exponential polynomial in 
f, of the same order as P(f), also with purely imaginary dominant 
characteristic roots. Let a(f) P(sin(f), cos(f)), r 2 rmfa > 
0 and ^(f) = e-'"^^l’^F{fmla) = e-’'^*P 2 (f). 


Multiplying the first identity through by m? and summing, we 
have 

Pnsin(nf -F p 2 ){n^ — m?) = 0. 

By the premise of the Lemma, B 7 ^ 0, so n(n — m) (n -F m) = 
0 , which is a contradiction. 

• Finally, only the symmetric cases Mi < M 2 = 0 and 0 = 
Ml < M 2 remain. Without loss of generality, by replacing / 
by — / if necessary, we need only consider the case 0 = Mi < 
Ms. 

Thus, assume 0 = Mi < M 2 . We now move our attention to 
the possible forms of P 2 . Since / has order at most 8 , it follows that 
P 2 has order at most 3. Thus, there are three possibilities for the set 
of dominant characteristic roots of P 2 : { 0 }, {±ic}, or { 0 , ±ic}, 
for some positive c G R IT A. We consider each of these cases in 
turn. 

First, if P 2 only has the real dominant eigenvalue 0, then P 2 is 
ultimately positive or ultimately negative, depending on the sign of 
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the most significant term of F 2 . Ultimate positivity of F 2 entails 
ultimate positivity of / as well, since P(sin(f), cos(f)) > 0 ev¬ 
erywhere, whereas an ultimately negative F 2 makes / change sign 
infinitely often. 

Second, assume the dominant characteristic roots of F 2 are 
{±ic}, so that 

f{t) = P(sin(f), cos(f)) + [D cos(cf -|- (^ 53 ) + 

for some ra > 0 and 1 /P 3 € M such that € A. Without loss 
of generality, we can assume c ^ Q, since otherwise, we are done 
by Theorem]^ But by LemmalU it will happen infinitely often that 
P(sin(f),cos(f)) = 0 and Dcos{ct -|- (fis) < —\D\/2, say. For 
large enough such t, < |P|/4, so we conclude that 

/ is infinitely often negative, and hence has infinitely many zeros. 

Third, assume the dominant characteristic roots of P 2 (f) are 
{ 0 , ±ic}, so that 

f{t) = P(sin(f), cos(f)) + e~’'^*(Pcos(cf -|- (^ 93 ) -|- E). 

We again assume c ^ Q, since otherwise the claim follows from 

Theorem[ 8 l Let M 3 E — |P| = mint>oP 2 (f). If M 3 > 0, 
then fit) clearly has no zeros. If M3 < 0, then there exists a 
non-trivial interval I C [0, 27r) such that if ct + ifis mod 27r G 
I, then P 2 (f) < 0. Since c ^ Q, Lemma |4] guarantees that 
P 2 (f) < 0 = P(sin(f), cos(f)) happens infinitely often, so / must 
have infinitely many zeros. Finally, if M 3 = 0, we argue that / 
is ultimately positive. Indeed, since P(sin(f),cos(f)) and P 2 (f) 
are both non-negative everywhere, /(f) = 0 can only happen if 
P(sin(f),cos(f)) = Dcos{ct + ipz) -f P = 0. This, however, 
would entail e** G A and G A, which contradicts the Gelfond- 
Schneider Theorem, since c ^ Q. Thus, we conclude / has no 
zeros. □ 

4.3 Three dominant oscillations 

Lemma 15. Let A, B, C, a, b, c be real algebraic numbers such 
that a,h,c> Q and A,B,C^ 0. Let also ipi,ip 2 ,>P 3 G R fie such 
that G A. Define the exponential polynomial f by 

f{t) = A cos(af + pi) + B cos(fef -I- P 2 ) + C cos(cf Fps) + D. 

It is decidable whether f has infinitely many zeros. 

Proof. The argument consists of three cases, depending on the 
linear dependencies over Q satisfied by a, b and c. 

Case I. First, if a, b, c are linearly independent over Q, then the 
claim follows directly from Theorem|7] 

Case II. Second, suppose that a, b, c are all rational multiples of 
one another: 

6 = — a, c = —a where n, m,k,l G N'*'. 
m I 

We make the change of variable t —>■ tml to obtain: 

fit) = Acosiiat)ml + pi) + B cos((af)ni -|- P 2 ) 

-I- C cos((af)fcm + ps) + D 
= P(sin(af), cos(af)), 

where P G A [a;, i/] is a polynomial obtained using the standard 
trigonometric identities. It is now clear that / is periodic, so it has 
either no zeros or infinitely many zeros. Let 

Ml min Pix,y) — min/(f), 

M 2 max Pix,y) = max fit). 

= l t>0 

Using the same reasoning as in Lemma [741 we see that Mi, M 2 
are algebraic and effectively computable: simply construct defining 


formulas in the first-order language C of real closed fields, and then 
perfo rm quantifier elimination using Renegar’s algorithm iRenegail 
Il992h . Then / clearly has infinitely many zeros if and only if 
All ^ 0 ^ AI 2 . 

Case III. Finally, suppose that a, b, c span a Q-vector space 
of dimension 2 , so that a, b, c satisfy a single linear dependence 
am + bn + cp — 0 where m,n,p G Z are coprime. At most 
one of the ratios a/ 6 , a/c and b/c is rational (otherwise we have 
dim span{a, b, c} = 1 ), so assume without loss of generality that 
a/c ^ Q and b/c ^ Q. 

Define the set 

T {x G [0, 27 r)® I Vm C 1 ? .u ■ ia, b, c) G 27 rZ => m ■ x G 27 rZ} 
= {(xi, X2, X3) G [0, 2 n)^ I mxi + nx2 + pxs G 27 rZ} 

Notice that if mxi -I- nx2 + pxz = 2 A: 7 r for some xi, X 2 , X 3 , then 
k < |m|-|-|n|-|-|p|,soT partitions naturally into finitely many 
subsets: T = U^i Tfe. where 

Tfc {(xi,X 2 ,X 3 ) G [0, 27r)® | mxi -|-nx 2 -l-pxs = 2fc7r} . 

def 

Consider the trajectory fi,(f) = {(af, bt, ct) mod 27r | f > 0}. De¬ 
fine also the sets R {fi,( 2 A; 7 r) | G N} and H {hit) \ t > 
0}. Because of the linear dependence satisfied by a, b, c, it is easy 
to see that P C P C T. By Kronecker’s Theorem, P is a dense 
subset of T, so clearly H must be a dense subset of T as well. 

Now define the function 

P(xi, X 2 , X 3 ) Acos(xi -I- pi) + B COs(x 2 -I- P 2 ) 

+ C cos(x3 F Pi) F D, 

so that the image of / is exactly {P(xi, X 2 , X 3 ) | (xi, X 2 , X 3 ) G 
H}. Let also the extrema of P over T be: 

Ml minP(xi, X 2 , X 3 ), 

T 

M 2 maxP(xi,X 2 ,X 3 ). 

Both of these values are algebraic and can be computed using 
quantifier elimination in the first-order language L of the real 
numbers: just use separate variables for cos(xj), sin(xj) and ap¬ 
ply the standard trigonometric identities to convert the linear de¬ 
pendence on xi,X 2 ,X 3 into a polynomial dependence between 
cos(xj), sin(xj). 

Now, by the density of H in T, if Mi < 0 < M 2 , then / must 
clearly be infinitely often positive and infinitely often negative, so 
it must have infinitely many zeros. The case Mi < 0 = M 2 
is symmetric to 0 = Mi < M 2 (just replace / and P by — / 
and —P, respectively), so without loss generality, we can assume 
0 = Ml < M 2 . In this case, we argue that / has no zeros, that is, 
even though P vanishes on some points in T, none of these points 
appear in the dense subset Ft. Indeed, consider the set 

Z {(cos(xi), sin(xi),..., cos(x 3 ), sin(x 3 )) | 

(xi, X 2 , X 3 ) G T, P(xi, X 2 , X 3 ) = 0} . 

Note that Z is clearly semi-algebraic, as one can directly write 
a defining formula in C from P(xi,X 2 ,X 3 ) = 0 and mxi + 
nx2 F PX3 G 27rZ. Moreove r, by the Zero-Dimensionality Lemma 
iOuaknine and Worrelll2014l Lemma 10), the function P(xi, X 2 , X 3 ) 
achieves its minimum Mi = 0 at only finitely many points in T*,, 
for each k. Since T is the union of finitely many T^, we imme¬ 
diately have that Z is finite. By the Tarski-Seidenberg Theorem, 
projecting Z to any fixed component will also give a finite, semi- 
algebraic subset of R, that is, a finite subset of A. Thus, we have 
shown that if P(xi, X 2 , X 3 ) = 0, then G A for all j — 1, 2, 3. 
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Now if f{t) = 0 for some f > 0, then we must have £ A, 

which by the Gelfond-Schneider Theorem entails ajc £ Q, a con¬ 
tradiction. □ 


4.4 One repeated oscillation 

Lemma 16. Let A, B, C, D, a, r be real algebraic numbers such 
that a,r > 0 and A ^ 0. Let also (pi,v52 £ M. be such that 
^ g®V 2 ^ ^ Define the exponential polynomial f by 

f{t) = t{Acos{at+(pi)+B) + {C cos{at+(p 2 ) + D)+e~'~^F{t) 

where F is an exponential polynomial with purely imaginary dom¬ 
inant characteristic roots. Suppose also that f has order at most 8 . 
It is decidable whether f has infinitely many zeros. 


On the other hand, we have the following chain of inequalities 
2 

> Sj } 

Ml! 

2 

= { definition of 5j } 

m\ + \D\) 

tj-i 

> { by f > tj-i } 

‘2{\C\+\D\) 

t 

> { triangle inequality and | cos (a;)! < 1 } 

ICI + |D| \Ccos{at + P2) + D\ 

t t ■ 


Proof. Since / has order no greater than 8 , it follows that F has 
order at most 2. Therefore, F{t) must be of the form E cos{bt-\-ip 3 ) 
for some i?, 6 £ R n A, 6 > 0, such that a/h ^ Q, and some 
such that £ A, since otherwise the imaginary parts of the 
characteristic roots of / are pairwise linearly dependent over Q, so 
our claim is proven immediately by Theorem[ 8 ] 

Consider first the magnitudes of A and B. If |j 4| > |B|, then 
the term tAcos{at + (pi) makes / change sign infinitely often, 
whereas if |B| > | A|, then for t large enough, the term tB makes / 
ultimately positive or ultimately negative, depending on the sign of 
B. Thus, we can assume I = |B|. Dividing / by B, and replacing 
pi by (pi +7r if necessary, we can assume the function has the form: 

f(t) = f(l — cos(af + pi)) 

{C cos{at P 2 ) + D) e~'"^E cos{bt ps). 

Considering the dominant term, it is clear that / is infinitely often 

def def 

positive. Let a{t) = t{l — cos(of -|- pi)), /3(f) = C cos(af -I- 

P2) + D and 7 (f) cos{bt -I- (ps). 

We now focus on the sign of the term /3(f) at the positive 

def 

critical times tj = —pi fa -|- 2jnfa (j £ Z) when 1 — cos(af + 

pi) vanishes. Notice that P{tj) = Ccos((p 2 — pi) + B 
M is independent of j. First, if M <0, then for all tj large 
enough, f{tj) < 0 , so the function must have infinitely many 
zeros. Second, if M = 0, then by the linear independence of a, b 
and Lemma |4l we have a{tj) = fi{tj) — 0 > 7 (fj) for infinitely 
many tj, so we can conclude / has infinitely many zeros. 

Finally, suppose M > 0. We will prove that / is ultimately 
positive. For each tj, define the critical region [tj — 5j,tj -1- (5j], 
given by 

^ def 2 ^|C| +|B| 

^ Oy/tj-l 

From here onwards, we only consider f large enough for any two 
adjacent critical regions to be disjoint. The argument consists of 
two parts: first we show /(f) > 0 for all large enough f outside all 
critical regions, and then we show /(f) > 0 for large enough f in a 
critical region. 

Suppose f is outside all critical regions and let j minimise 
|f — tj\. Since the distance between critical points is 2it/a by 
construction, we have a|f — f^j < it. Therefore, 


Combining, we have 
a(f) + ^(f) > a(f) - |^(f)| 

= f(l — cos(af + pi)) — |C'cos(af + P 2 ) -b D\ 
>\C\ + \D\. 

Thus, if f is large enough to ensure | 7 (f)| < |C| + |B|, we have 
/(f) > 0 outside critical regions. 

For the second part of the argument, we consider f in critical 
regions. Notice that the values of /3(f) on \tj — Sj,tj -|- are 
independent of the choice of tj. Moreover, we have j3{tj) = M > 
0, so there exists some e > 0 such that for all f £ [tj — e, tj -|- e], we 
have /3(f) > M/2, say. Now for any critical point tj chosen large 
enough, we will have [tj — 5j, tj -b ] C [tj — e,tj + e], so (f) > 
M/2 on the entire critical region. Let also tj be large enough so 
that for any f in the critical region, we have | 7 (f)| < M/2. Then 
wehave/(f) = a(f)-b/3(f)-b7(f) > /9(f) —| 7 (f)| > 0 ,completing 
the claim. □ 

5. Hardness at Order 9 

Diophantine approximation is a branch of number theory con¬ 
cerned with approximating real numbers by rationals. A central role 
is played in this theory by the notion of continued fraction expan¬ 
sion, which allows to compute a sequence of rational approxima¬ 
tions to a given real number that is optimal in a certain well-defined 
sense. For our purposes it suffices to note that the behaviour of the 
continued fraction expansion of a real number a is closely related to 
the following two constants associated with a. The Lagrange con¬ 
stant (or homogeneous Diophantine approximation constant) of a 
is defined by 

{ \ Tl \ C 'I 

c : a-< —IT for infinitely many m, n £ Z > 

I m I m-^ J 

By definition Loo (a) is a non-negative real number. 

A real number a is called badly approximable if Loo (a) > 
0. The badly approximable numbers are precisely those whose 
continued fraction expansions have bounded partial quotients. 

Khinchin showed in 1926 that almost all real numbers (in the 
measure-theoretic sense) have Lagrange constant equal to zero. 
However, information on the Lagrange constants of specific num¬ 
bers or classes of numbers has prov en to be el usive. In particular, 
concerning algebraic numbers, Guy jGu ^l 200 l) asks 


2 


< 1 — cos(af — atj) = 1 — cos(af -\- pf). 


Is there an algebraic number of degree greater than two 
whose simple continued fraction expansion has unbounded 
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partial quotients? Does every such number have unbounded 
partial quotients? 

The above question can equivalently be formulated in terms of 
whether any algebraic number of degree greater than two has stricly 
positive Lagrange constant or whether all such numbers have La¬ 
grange constant 0. 

Recall that a real number a is computable if there is an algorithm 
which, given any rational e > 0 as input, returns a rational q such 
that |q — a:| < e. We can now state the main result of the section. 

In this section, we will show that a decision procedure for the 
Infinite Zeros Problem would yield the computability of Lao {a) for 
all a £ R n A. 

Fix positive a£RnA, c£Q and define the functions: 
fi{t) e*(l — cos(?)) — cos(a?)) — csin(a?), 

f 2 {t) e*(l — cos(?)) -F t{l — cos(a?)) + csm{at), 

f{t) e*(l — cos(?)) + t{l — cos{at)) — c\ sin(o?)| 

= min{/i(f),/ 2 (?)}. 

It is easy to see that /i and /2 are exponential polynomials of order 
9, with six characteristic roots: three simple (1 and 1 ± i) and three 
repeated (0 and Lai). Thus, the problem of determining whether 
fj has infinitely many zeros is an instance of the Infinite Zeros 
Problem. Moreover, it is easy to check that / has infinitely many 
zeros if and only if at least one of fi and /2 has infinitely many 
zeros. 

We will first state two lemmas which show a connection be¬ 
tween the existence of infinitely many zeros of / and the Lagrange 
constant of a. We defer the proofs to Appendix I aI 

Lemma 17. Fix a £ R n A and e,c € Q with a, c > 0 
and e £ (0,1). If f{t) = 0 for infinitely many t > 0, then 
Lao{a) < c/27r^(l — e). 

Lemma 18. Fix a £ R n A and e, c £ Q with a, c > 0 and 
e £ (0,1). If Loo{a) < c(l — e)/ 27 r^, then f(f) = 0 for infinitely 
many t > 0. 

We now use the above lemmas to derive an algorithm to com¬ 
pute Loo (a) using an oracle for the Infinite Zeros Problem, estab¬ 
lishing our central hardness result: 

Theorem 19. Fix a positive real algebraic number a. If the Infinite 
Zeros Problem is decidable for instances of order 9, then Lao (a) 
may be computed to within arbitrary precision. 
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/(?) e*(l — cos(i)) + t{l — cos(at)) — c| sin(a?)| 
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Lemma 17. Fix a £ R n A and £,c € Q with a, c > 0 
and £ £ (0, 1). /f f{t) = 0 for infinitely many t > Q, then 
Aoo(a) < c/27r^(l - e). 


Proof. Suppose f{f) = 0 for infinitely many t. Clearly, this also 
entails f{t) = 0 for infinitely many t > T, for any particular 
threshold T > 0. (Indeed, f{t) = min{/i(f),/ 2 (f)} for expo¬ 
nential polynomials /i and /2 given at the beginning of Section|3 
Thus, on any bounded interval, / has no more zeros than /i and /2 
combined, i.e., only finitely many, by the analiticity of /i and / 2 .) 
We will show that T can be chosen in such a way that every zero 
of f{t) on [T, 00 ) yields a pair (n, m) £ which satisfies the 
inequality 

Ini c 

I m I ^ 27r^m^(l —e) 

This is sufficient, since infinitely many zeros of / yield infinitely 
many solutions, and therefore witness Loo(a) < c/27r^(l — e). 

Thus, consider some t such that f{t) = 0 and t > T for 
some threshold T to be specified later. Let t = 2nm -|- and 
at = 27rn -I- S 2 , where m, n £ N and di, ^2 € [—it, it). Then we 
have 

I n I _ 1^2 — ci5i| 

I m I 2itm 

We will show that for T large enough, f{t) = 0 for f > T allows 
us to bound | S 2 \ and | a5i | separately from above and then apply the 
triangle inequality to bound 1^2 — a5i|. 

First, choose ipi, (^2 G (0,1) such that 1 —952 > 1 —(pi > 1—e. 
Let T be large enough for the following property to hold: 


t + 1t ^ 1 — ifi2 
t — 2it ~ 1 — tpi 


for all t > T. 


In particular, since m = {t — Si)/2it and |5i | < tt, we have 


2m ^ t + it ^ l-ip2 
2m —1 ~ t — 2it ~ 1 — ifii 

Let also T be large enough to make the following property valid: 

1 — cos(®) < A |x| < TT ^ < 1 — cos(a;). (6) 

Now we have the following chain of inequalities: 


1 — cos(52) 

< { f{f) = 0, noting e*(l — cos(f)) > 0 } 
c| sin((? 2 )| 

t 

< { by I sin( 2 :)| < \x\ } 

c\S2\ 

t 


Then by llbll, we have 


1 — cos((52) > 


(1 - ip2)Sl 


Thus, combining the upper and lower bounds on 1 — 005 (^ 2 ) and 
using a on the last step, we have 


\52\ < 


2c 


< 


2c 


< 


t{l — 1 P 2 ) {2m — l)it{l — LP 2 ) mit{l — Lpi) 


Second, let q: (1 — e) —(1 —:^i) ^ > 0. Let the threshold 


T be large enough so that 


-t , coi 
e < 


27r 


dTT^O^ \t + It 


for t > T 


(7) 


and 


if 1 — cos(a;) < c/e^ and |a:| < it, then a:^/4 < 1 — cos( 2 ;). 


The following chain of inequalities holds: 
1 — cos(5i) 

= {by/(f) = 0} 
c| sin((52)| — t{l — cos((52)) 


(8) 


< { by I sin(52)|, | cos(52)| < 1} 


< { by 0 } 

/ 27r 

4:it^a^ y f -|- TT 

< { by |5i| < TT } 
27r 


— (5i _ 

= {t = 2itm + } 


Moreover, as 1 — cos(5i) < ce“* < ce~^, by 0, we have 
1 - cos(Oi) > —, 

so combining the lower and upper bound on 1 — cos((5i), we can 
conclude 

|a5i| < —. 
itm 

Finally, by the triangle inequality and the bounds on |a(5i| and 
|(52|, we have 




1^2 — a5i| ^ 1^21 + |a5i| 


< 


2itm 

c 


a + 


2itm 
1 


1 — (fi J 27r^m^(l — e) 


Now, by the premise of the Lemma, there are infinitely many t > T 
such that f{t) = 0, each yielding a pair {n,m) £ which 
satisfies the above inequality. These infinitely many pairs (n, m) 
witness Lao{a) < c/2it^{l — e), as required. □ 

Lemma 18. Fix a £ R n A and e, c £ Q with a,c > 0 and 
£ £ (0, 1). IfLoo{a) < c(l — e)/27r^, then f{t) = Q for infinitely 
many t. 


Proof. We will show that there exists an effective threshold M, 
dependent on a, c, £, such that if 



< c(l -e) 

~ 27r^m2 


(9) 


for natural numbers n, m with m > M, then f{2itm) < 0. Note 
that this is sufficient to prove the Lemma: the premise guarantees 
infinitely many solutions (n, m) £ of 0, so there must be 
infinitely many solutions with m > M, each yielding f{2itm) < 
0. Since / is continuous and moreover is positive for arbitrarily 
large times, it must have infinitely many zeros on [27rM, oo). 

Now let M be large enough, so that c(l — £)fitM < it and 


if I®! < c(l — £)/itM, then (1 — e)|x| < | sin(a:)|. (10) 

Suppose that 0 holds for n, m £ N with m > M and write 
t = 2itm. We will show that f{f) < 0. By 0, we have \am — 
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n-l < c(l — E)l2-K^ra. Therefore, at = 2'Kam = 27rn + S where 
1^1 < c(l — e)l'Km < TT. We have 

/(i) 

= { as cos(f) = 1 } 
t{l — cos(5)) — c| sin(5)| 

< { by i fTOt and 1 — cos{x) < ®^/2 } 

7rm5^ — c(l — e)|(5| 

< { by |5| < c(l — £)/7rm } 

0 . 

□ 
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